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Abstract. It is known, due to Mordukhai-Boltovski, Ritt, Prelle, Singer, Christopher and 
others, that if a given rational ODE has a Liouvillian first integral then the corresponding 
integrating factor of the ODE must be of a very special form of a product of powers and 
exponents of irreducible polynomials. These results lead to a partial algorithm for finding 
Liouvillian first integrals. However, there are two main complications on the way to ob- 
taining polynomials in the integrating factor form. First of all, one has to find an upper 
bound for the degrees of the polynomials in the product above, an unsolved problem, and 
■ then the set of coefficients for each of the polynomials by the computationally-intensive 

method of undetermined parameters. As a result, this approach was implemented in CAS 
only for first and relatively simple second order ODEs. We propose an algebraic method 
| for finding polynomials of the integrating factors for rational ODEs of any order, based on 

OO ' examination of the resultants of the polynomials in the numerator and the denominator of 

the right-hand side of such equation. If both the numerator and the denominator of the 
right-hand side of such ODE are not constants, the method can determine in finite terms 
an explicit expression of an integrating factor if the ODE permits integrating factors of the 
above mentioned form and then the Liouvillian first integral. The tests of this procedure 
based on the proposed method, implemented in Maple in the case of rational integrating 
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■ factors, confirm the consistence and efficiency of the method. 



Key words: differential equations; exact solution; first integral; integrating factor 
2000 Mathematics Subject Classification: 34A05; 34A34; 34A35 



X 

1 The Prelle— Singer method 

u is an integrating factor for an nth order ODE in solved form 

-r^ - f(x,yo,yi,--- ,y n -i) = 0, (l) 

where yo = y(x), Vj = ^-j^r > by standard definition pQ, if /j,(y n — f) is a total derivative of some 
function ((x, y , yi, . . . , y n -i), that is, 



/' 



d n y ,\ d( U d( 

-d^- f ) = d- X + ^ +1 dy- J - (2) 



Let 



V = a • ( 3 ) 



dy. 



n-l 
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Then © becomes 
■ c 



^(C) + /tt^ = 0, (4) 



where 

D = — + V — 

We conclude that the integrating factor of ODE ([T]) is = g^C _ ; where the function £ is a first 
integral, a solution of the linear first-order PDE © . 

If we eliminate the function £ from the system (J3J), (@J), we obtain a PDE system for the 
integrating factor. One of the equations of the PDE system obtained by double differentiation 
of (J3J) with respect to y n -\ is as follows (n > 2): 



fy* \ + d 2 (i» + o dfx 
dy n -ij dy 2 n -\ dy n 



Assuming that an integrating factor is known, it is well-known that we are able to obtain cor- 
responding first integral by quadratures. If we know n independent first integrals for ODE (^Q), 
we can then find its general solution, at least in implicit form. As it is seen from ©, the alge- 
braic structure of the integrating factor is simpler than the structure of the first integral, so in 
some cases finding the integrating factors is an easier problem. 

The most beautiful way of finding integrating factors is the Darboux method and its refine- 
ments. If / in (^Q) is a rational function with respect to all variables x,yo,y\, . . . ,y n -i, we will 
call such an ODE a rational one. In 1878, Darboux |2j proved that if a rational ODE has at 
least m(m+ l)/2, where m = max(order (numerator (/), order(denominator(/)) invariant algeb- 
raic curves Pj, now known as Darboux polynomials, then it has a first integral or an integrating 
factor of the form ]J Li P" -1 for suitable complex constants a,. So, in principle, the problem of 
finding an integrating factor is reduced here to the problem of finding Darboux polynomials. 

In 1983, Prelle and Singer [3] proved the essential theorem that all first order ODEs, which 
possess elementary first integrals, have an integrating factor of the above mentioned form JX P-°*. 
Their result unifies and generalizes a number of results originally due to Mordukhai-Boltovski [1] , 
Ritt jS] and others. Singer in jHj refines this result for ODEs, which possess Liouvillian first 
integral. 

The following theorem is a corollary of the Proposition 2.2 |Sj and of the definition of the 
integrating factor © (see also [Jj): 



Theorem 1. If the n-th order ODE 

d n yp = A(x,y ,y 1 , . . .,y n -i) 
dx n B(x,y ,y 1 ,...,y n - 1 y 



(6) 



where A,B £ K[x,yo,y\, . . . ,y n -i] are polynomials, has a local Liouvillian first integral, then 
there exists an integrating factor of ODE JBJ) of the form 



p = Y[P?exp[b l[C#), (7) 



where Pi 's, Qj 's 6 K[x,yo,yi, . . . ,y n ~i] ar & irreducible polynomials, ai 's, &o are constants and 
bj 's are integers. 
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So, to find an integrating factor to ODE (JHJ) in the conditions of Theorem ^ we first should 
obtain the sets of polynomials Pi and Qj, and if we succeed in this task then we can obtain the 
constants cjj's and bj's from the full PDE system for the integrating factor. 

To consider the ways to find the sets of polynomials Pi and Qj, let us substitute the integrating 
factor in the form (J7J) into (|5]>. After some rearrangement, we arrive at 



UB = 2A 



dB 

dy n ~ 



n*n«: 



where L\ is a polynomial. As A and B are relatively prime polynomials, we conclude that 
B divides one of the polynomials P, or Qj. Suppose that 

A = Ai(y n _i)A 2 (x,y ,yi,. . . ,y„_i), 

B = B 1 (y n - 1 )B 2 (x,y ,yi, . . . ,y n - 2 )B 3 (x,y ,yi, ■ ■ ■ ,y n -i) 

and 



Substituting these expressions for A, B and \x into (|Kj). we can conclude, by considering divis- 
ibility of polynomials, that it is necessary that £3 = 1 and 773 = 0. That is, if - = or 

D{Pi) = ( - = or D(Qj) = 0) then such polynomials are factors of A or B, and they can 
be easily selected and included into the set of candidates for the integrating factor structure. 
Without loss of generality, we will consider an integrating factor of the form 



l-L = B\\P^^ v [b \[Q) 



(8) 



and will focus our attention on the cases when J^ Pi - ^ and D(P{) 7^ ( q®' 1 7^ and 
D{Qj) + 0). 

For the case of 7^ 1, substitution of (jSJ into © leads to 



L 2 p k = a k (a k - i)^—*- n ^ 2 n ^ 2 



BD(f5fc) + A 



dP k 
dy n - 



(9) 



and similarly for the case of b k < 



n-l 



ip?n«. 



2(6j+l) 
J 





3^k 



BD(Q k )+A 



dQk 
dy n -\ 



(10) 



where L 2 and L3 are polynomials. 

All components involved in © and (|1(J|) are polynomials. We assume that Pi and Qj are 
irreducible, relatively prime polynomials, so gcd(Pt ; Qy Pi - ) and gcd(Qj, g^' t ) are constants, 
where greatest common divisor is denoted by gcd. We can conclude that (/ | g means that / 
divides g) 



Pi I BD(Pi) + A 



OP 



Oy 



n-l 



(11) 
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Qj\BD{Qj) + A^-. (12) 
oy n -i 



and 



The basic significance of relations 1)1 lj) and (|12|> is that we can investigate each of the poly- 
nomials Pi or independently. Properties (|11|) and 1)12|) constitute the basis of the so-called 
Prelle— Singer method, which originates from Darboux [2], to obtain the integrating factors and 
first integrals of ODEs of type ©• Supposing that 

N N N 

ft=EE'" E ^O-in-!^ -^-! 1 (13) 
j=0 jo=0 Jn-1=0 

and substituting into (|lljl. we can in principle obtain Ckij -j n _ 1 , s and as a result Pfc (or 
similarly Q&) in explicit form. 

To do so, we first have to establish an upper bound N for the degrees of the polynomial Pi 
(or Qk)- It is known that ./V exists, but there is no constructive method to evaluate it so far. It 
is the principal weakness of the Prelle-Singer method. In the existing implementations of the 
Prelle-Singer method, an upper bound N is predefined by a user of the method. 

With an increase of equation order the number of undetermined variables Ckijo---j n 

, and, 

correspondingly, the dimension of the system of algebraic equations for these constants increases. 
As a result, solving of this system becomes problematical. This approach was implemented in 
CAS only for first and relatively simple second order ODEs |8] \9\ HHj ITT ) IT2l IT3] HI) ITK] . and as 
a rule the method can not be used for N > 4. In contrast, the modified method proposed in 
the next section does not suffer from these restrictions. 

We must also mention that a rational ODE does not always have a local Liouvillian first 
integral. While generalizations of the Darboux integrability theory allow in some ( "integrable" ) 
cases to find non-Liouvillian first integrals (see, e.g. El El El ED]), we consider only 
Liouvillian first integral theory. 

2 The method of resultants 

Let us investigate closely the relation ()TT]) or equivalently (JT2J). Let R z (f,g) be the resultant of 
polynomials / and g with respect to the indeterminate z. In the sequel, we will consider only 
the case when both A and B are not constants and there are such z € (x, yo, y\, . . . , y n -i) that 
R Z (A, B) is not a constant. 

The following Theorem [2] and Corollary are our main results: 

Theorem 2. If a polynomial Pi satisfies the condition (|1 1ft then for any indeterminate z € 
(x, yo, yi, . . . , y n -i) the resultant R z (Pi, B) must divide the resultant R Z (A, B). 

Proof. We may write Pi L = BD(Pi) + A Q® Pi - for some polynomial L. Taking resultants of 
both sides of this equation with B and with respect to an indeterminate z 6 (x,yo,yi, ■ ■ ■ , y n -i), 
we get the following polynomial equation 

R z {Pi,B)R z (L,B) = R Z (A,B)R Z (^-,B 



dy n - 

As Pi and Qy ^ have no common roots, R z (Pi,B) does not divide Rz\ Q y ' 1 ; Bj . Therefore, 
R z (P l ,B)\R z (A,B) (14) 

and it is obvious that we can choose R z (Pi,B) as a divisor of R Z (A,B) up to a constant fac- 
tor. ■ 
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How can we recover Pj if we know the resultant R z (Pi, B)l First of all, we note that the 
case when A/B depends on only one indeterminate is trivial. As we have mentioned above, if Pj 
depends only on one indeterminate, then such polynomials are factors of A or B and we have 
already included them into the set of candidates. So, we assume that Pi depends at least on two 
indeterminates and there exists such z that R z (Pi,B) is a polynomial. We do not insist that 
R z (Pi,B) is irreducible. 

It is known from elementary facts about resultants that there exist such polynomials aj,/3j £ 
K[x,yo,yi,...,y n -i] that /3jPj = o-j B + R z {P, h B). 

As a corollary of Theorem[21 we have the following result. Assuming that Pj's (and Qj's) are 
irreducible polynomials we conclude that 

Corollary. If polynomials Pi (and respectively Qj) satisfy (|14j) then the following holds 

P« hyp = onB + a [one of nonconstant factors of R Z (A, B)], (15) 

where Pihyp = ftiPi, A is some polynomial, cti € ( — 1,0,1) and Cj are constants. 

Proof. If ai ytz identically, then for any z € (x, yo, y\,..., y n -i) from 

R z ((3iP l -a i B,P i )= 1 R z (a i ,P l )R z (B,P l ) 

= CjP 2 ([one of nonconstant factors of R Z {A, B)],B) = R z (B,Pi), 

where 7 G (—1,1), we conclude that R z (a.i,Pi) = I/7. So as Pj depends at least on two 
indeterminates then «j must be a constant and Oj £ (— 1,0,1). 

Now we are able to modify the Prelle-Singer method by replacing the hypothesis l)13J) 
with ((To]) and finding constants en and Cj from the following equation 

R z ^BD(P ihyp ) + A^^,P ihy ^j = 0. 

Then Pj is a irreducible factor of Pjhyp which satisfies (|11[) . 

Thus, we can use algebraic, not differential, relations for obtaining polynomials as in the 
original Darboux (and Prelle-Singer) approach. First of all, we bypass the principal weakness 
of the Prelle-Singer method as we do not need to establish an upper bound N for the degrees 
of polynomial Pj here. Second, we have the advantage that finding only two undetermined 
constants a« and Cj in the Pjhyp structure is much easier than finding the set of Ckij —j n _ 1 , s. ■ 

Example. Let us consider a relatively simple second order ODE 

„ = y' 3 + y'(x - 2) - y 
y y v + (2y> -l)(x + y )- x - 

For this ODE, R X (A,B) = —(2yf + y\ — 2)(yf — 2y\ — yo) and we can form the following 
hypothesis Ph yp = a(yf + (2yi — l)(x + yo) — x) + c{y\ — 2y\ — yo). Calculation of the following 
resultant leads to R yi (BD(Ph yp ) + A ^g^ p , Phyp) — a 2 (a — c)L = 0, where L is some polynomial. 
So, there are two cases, when a = then Pihyp = c(y 2 — 2y\ — yo), and when a = c then 
^2 hyp = 2c(yi - l)(yi + y + x). 

By substitution of the irreducible factors of Pi h yP and P2 h yP to (fTTj) , we select the following 
candidates for the integrating factor structure: Pi = (yf — 2y\ — yo) and P2 = {y\ + yo + x). So 
the hypothesis of /x is 



Mhyp = (Vi + (2yi -l)(x + y ) - x)(y 1 + y + x) xi (y\ - 2y 1 - y ) 
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After substitution of this //hyp into the PDE system for the integrating factor, we obtain that 
XI = —(XI + 2) so an integrating factor of given ODE is 

fi = (y> 2 + (2y> - l)(x + y) - x)(y' + y + x) xi (y' 2 - 2y' - y)^ xl+2 \ 

where XI is an arbitrary constant. 

Examinations of the cases when a k = 1 and b k > 0, are more difficult. We are not able to 
give here the full analysis of the problem which necessitates consideration of many subcases. 
For brevity, we will only demonstrate that the method of resultants in principle enables us to 
obtain all the needed polynomials in a finite number of steps. The main tool here is usage of 
repeated resultants. 

For compactness, we will use the following notation for the repeated resultants 



,(/, (g,h)) =R Z2 (R Z1 (f,h),R Zl (g, h)). 



For the case = 1, we obtain the following equation similar to @: 

up k = n ^ n qj( m i a + m ? b ) 



i^k 



* i 
i^k 



2^- + D(B) 

oy n -i 



m + n(P k y SB 



dy n 



dy n - 



where L4, Mi, M2 are some polynomials, and Pj (i ^ k), Qj are only such polynomials from 
the integrating factor structure which do not satisfy conditions and (|12|). So for this case 
we can obtain the following properties ( Q y Pk - 7^ 0): 



So«(ft,WB)) I So, 



( II' II. 

\ i j 
i+k 

+ D{P k 



dy n - 



2^-+D(B) 

OVn-l 

(ABU 



dP k 



dy 



n-l 



and further 
R 



21,22,23 



Pk,[A,B, 



OB 

dy n - 



R 



21,22,23 



II' II'- 

. i 
i^k 



dP k 
3 dy n -i 



; 8A 

' dy n - 



+ D{B) 



A,B, 



dB 

dy n ~ 



(16) 



(17) 



From (|16|) we can obtain the hypothesis for P k and then with the help of Q17|) all its unde- 
termined constants. 

If there is a polynomial P s in the integrating factor structure which satisfy (jllj) . then 
finding P k can be simplified. Since 



R z (P k ,P s )\R z [A\\pf\{Q) 

\ i j 
i^k 



dP s 

dy n - 



dP 

Aj-^ + BD{P k ) 
dy n -i 



,P.\, 



(18) 
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then 



R Zl , Z2 (Pk, (Ps,B)) | R ZUZ2 U 2 H Pf H Q) 

\ i j 




2 dP s dP k 



(Ps,B) 



) 



(19) 



Among other things we can observe from ()16|) and (|18|) is that it is very likely that the P^'s 
satisfy the relation (jllj) even if = 1. Also, there are relations for Qk with b)~ > 0, similar to 



The main elements of the method described above were implemented in demonstration pro- 
cedures in Maple [211 I22| for the case of rational integrating factors. 

We have tested routines in the following area: for a given ODE order we assigned the first 
integral £ of the following type 



where Ni, Mi, ^ are some polynomials, oti > 0, > 0, r/i, 9i are some constants. Such 
algebraic first integral leads to an ODE of type Q with / = —D{Q)/ ^ k - , which is potentially 

tractable by the proposed procedure as this ODE has a rational integrating factor fi = = 



The first stage - finding candidates and hypothesis of the structure of an integrating factor - is 
very fast. The second stage is relatively resource consuming - finding the unknown parameters d{ 
of the structure of an integrating factor. 

We can conclude that our prototypes are workable for rational ODEs even with symbolic 
constant parameters of orders from n = 1 to ?i = 4 — 5. Although in principle the method 
provides ability to find an integrating factor for almost any of the above mentioned type of 
ODEs, in practice the procedure heavily relies on some basic CAS functions such as factor, 
simplify and so on, and in some cases fails. In addition, different outputs are observed in 
different Maple versions for the same ODE. 

As a rule, Maple cannot find any integrating factor for ODEs from the considered test 
area while proposed procedure is usually successful. Often our procedure produces more than 
one integrating factor. For example, the output may contain a set of /x's or /i with arbitrary 
parameters, and sometimes they lead to independent first integrals. 

3 Conclusion 

We have proposed an algebraic method for finding polynomials Pi and Qj of the integrating fac- 
tor form (J?! of rational ODE © for any order using resultants of the numerator and denominator 
of the ODE's right-hand-side. 

The main advantage of the method of resultants is ability to obtain the polynomials Pi and Qj 
without an a priori assumption of the upper bound for their respective degrees. Our method 
bypasses the principal weakness of the Prelle-Singer (Darboux) method. 

Moreover, the number of undetermined variables involved in calculations by our method is 
radically smaller than in the Prelle-Singer (Darboux) method. This property allows us to obtain 
polynomials Pi and Qj of arbitrary degrees without any complications whereas the Prelle-Singer 
(Darboux) method can not be useful as a rule for N > 4. 
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